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A renormalization group treatment of metamagnetic quantum criticality in metals is presented. 
In clean systems the universality class is found to be of the overdamped, conserving (dynamical ex- 
ponent z — 3) Ising type. Detailed results are obtained for the field and temperature dependence of 
physical quantities including the differential susceptibility, resistivity and specific heat near the tran- 
sition. An application of the theory is made to Sr3Ru207, which appears to exhibit a metamagnetic 
critical end-point at a very low temperature and a field of order 5 — 7T. 



A metamagnetic transition is empirically defined as a 
rapid increase in magnetization at a particular value of 
applied magnetic field. Because there is no broken sym- 
metry involved, one expects a first order transition from 
a low magnetization to a high magnetization state as 
an applied magnetic field H is swept through a (tem- 
perature dependent) critical value H mm (T), The curve 
of first order transitions H mm (T) terminates in a crit- 
ical point (H*,T*). By appropriately tuning material 
parameters it is possible to reduce T* to 0, yielding a 
quantum- critical end-point. This situation is depicted in 
Fig. 0: (b) shows a typical metamagnetic line and critical 
end-point in the field-temperature plane and (a) shows a 
possible variation of the temperature of the critical end- 
point with pressure. 

A number of 'strongly correlated metals', including 
UPt 3 Q, CeRh 2 Si 2 g and other heavy fer mion com- 
pounds [ |] , as well as d-electron systems such as MnSi 
exhibit metamagnetic transitions with properties sugges- 
tive of proximity to a quantum critical point. Recent 
measurements strongly indicate that the bilayer ruthen- 
ate d-electron system Sr 3 Ru 2 07 [|| is, at ambient pres- 
sure and moderate applied field, tuned almost exactly to 
such a quantum-critical end-point ||. 

Metamagnetism in metals has not been extensively 
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FIG. 1: (a) Schematic phase diagram, showing variation of 
end-point of line of metamagnetic first order phase transitions 
as control parameter (e.g pressure) is varied, (b) Schematic 
phase diagram in H, T plane for p < po showing metamagnetic 
line and location of end-point. 



studied. Mean-field treatments have recently ap- 
peared ^] and some discussion was given in the con- 
text of a treatment of weak ferromagnets via the 'SCR' 
method ||, but the critical phenomena have not so far 
been investigated. Apart from the experimental rele- 
vance, the issue is of fundamental importance. In a metal 
at T — a nonzero magnetization corresponds to a split- 
ting of the "spin up" and "spin down" Fermi surfaces 
and a first-order metamagnetic transition corresponds to 
a jump in this splitting. At the metamagnetic quantum- 
critical end-point fluctuations of the magnetization, and 
therefore of the Fermi surface positions, become large: in 
a certain sense the material does not have a well defined 
Fermi surface and may therefore be thought of as a non- 
Fermi-liquid. In this paper we give what we believe is 
the first renormalization group theory of metamagnetic 
quantum criticality. We present general results, which 
may be applicable to a number of systems, and a detailed 
application to S^R^Oy. We use the standard approach 
to metallic quantum criticality JToj , O], which involves 
integrating out the electron degrees of freedom to ob- 
tain a model of overdamped bosonic excitations which 
are analysed by renormalization group methods. 

Quotes are placed about "spin up" and "spin down" 
above because in many metamagnetic materials spin or- 
bit coupling is large and spin is not a good quantum 
number. However, for most purposes one may adopt 
a 'pseudo-spin' notation [|l2| labeling the two Kramers- 
degenerate states in zero-field. The Kramers degeneracy 
is broken by an applied field, leading to two Fermi sur- 
faces and the theory carries through as in the non-spin- 
orbit case with one important exception noted below: if 
spin-orbit coupling is strong, impurity scattering may af- 
fect the dynamics differently. 

To treat the critical behavior we assume that the fluc- 
tuations in the relative Fermi surface positions may be 
represented by fluctuations of a magnetization variable 
to. We write an action for longitudinal fluctuations 
tP(x,t) — (|to| — m av (_ff*)) /toq of the magnetization 



2 



density m about its average value m av at the critical 
field H* , normalized to some typical magnetization den- 
sity mo (for example, the high-field saturation magneti- 
zation). We define the critical field H* by the require- 
ment that at T — the action has no static third order 
terms. We write h = (H — H*) /H*, introduce a cutoff 
length a (for example, the lattice constant) and define an 
energy scale E c by the requirement that the coefficient 
of the static quartic term is 1/4; the action in d space 
dimensions and imaginary time becomes 
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Here 5 (which may be varied e.g. by changing pres- 
sure) tunes the system through the metamagnetic critical 
point and Sd yn (discussed below) expresses the order- 
parameter dynamics. 

We have assumed in writing the static part of 5 me t a 
that the gradient expansion is the conventional one, that 
the coefficients are simple numbers and that the parame- 
ters vary with temperature only as T 2 , as usual in Fermi- 
liquid theory. However, several recent papers have called 
these assumptions into question |l4|, 15[ . In particu- 
lar, Bclitz, Kirkpatrick and Vojta \14\ have presented 
perturbation-theory results indicating that the expansion 
of the static spin susceptibility of a conventional, non- 
critical clean Fermi-liquid about the q = value contains 
terms of order \q\ (in d = 2). Such terms would in- 
validate the usual gradient expansion. Subsequently, a 
closely related \T\ temperature dependence of the sus- 
ceptibility and other parameters was studied [fl5|| . The 
available perturbative calculations |l4| suggest that spin 
rotational invariance is required to obtain the anomalous 
momentum dependence, so that in the present situation 
the gradient expansion might be the conventional one. 
However, one might expect £q to be unusually large if 
m av (H*) is small. A detailed investigation of the behav- 
ior near a metamagnetic point could therefore be a useful 
test of this still unsettled issue. 

The dynamic part S'dyn follows because the order pa- 
rameter is essentially the difference in position of the 
spin-up and spin-down Fermi surfaces. Fluctuations at 
nonzero q correspond to locally increasing the number 
of spin-up electrons and decreasing the number of spin- 
down electrons. If spin is conserved such a fluctua- 
tion can relax only via propagation or diffusion of elec- 
trons within each spin manifold. In a clean spin-orbit- 
coupled system, pseudo-spin is conserved (at least for 
fields aligned along a crystal symmetry axis) and the 
same arguments apply. Therefore, in a clean system 
one expects (the term is most conveniently written in 
frequency-momentum space) 
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corresponding to overdamped but conserved fluctuations, 
and yielding the dynamical exponent z — 3. Here v is a 
velocity, presumably of the order of the Fermi velocity 
and ip{q,w n ) = J ^E c dTe^ s - w " T ip(x,T). Note that 
we are concerned only with longitudinal fluctuations, 
so 'precession' terms d T ip ■ V 2 ip x ip are not important. 
Strong (pseudo)-spin-conserving scattering would lead to 
diffusion (\uj n \ jvq — > |a>„| / Dq 2 ) changing z to 4. A mo- 
mentum non-conserving spin orbit coupling (as from im- 
purities in the presence of strong spin orbit scattering) 
would lead to relaxation {i.e. Eq. || with vq replaced by a 
momentum independent scattering rate) implying z = 2. 
An important scale is the characteristic energy tu s { of a 
spin fluctuation at momentum q c = 1/a, u) s { = v^ 2 q 3 . 
Typical uj s f values for transition metal magnets are of 
the order of 500K; for heavy fermion systems they are at 
least an order of magnitude smaller Q . 

We analyse the theory by the usual one loop renor- 
malization group equations Jil"| which, after mode elimi- 
nation and rescaling, relates the theory with parameters 
5,u, h to a new theory with parameters 5',u',h'. The 
behavior at h = has been previously reported [ Q |; we 
focus here on the h dependence. The scaling equations 
are (we assume henceforth that z — 3) 
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The field h scales as h(X) = e^ A and T(A) = Te 3X . The 
effect of eliminated modes on 6 is contained in / which 
is calculated by expanding the theory about the value 
ip{5h) which extremizes the static part of S'mcta at the 
rescaled field, and then using the Gaussian approxima- 
tion to the resulting action to evaluate the integral over 
eliminated modes. Operationally, this means that we cal- 
culate / assuming the scales of interest are larger than the 
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running 'mass' r c ff = S + 3uip and then stop scaling at 
r c ff = 1 . Expressing momenta and frequencies in units of 
q c = 1/a and u st then / = A / coth(f) r ^ 



{2-K) d v \q\ 
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with A = (lu s [/E c ) (a/£o) 2 (the ' denotes summation over 
eliminated modes). 

The solution of Eqs. || and |^ follows and is dis- 
cussed in detail elsewhere [[l6|. Because in all cases of 
physical interest the model is at or above its upper crit- 
ical dimension, quantal fluctuations lead only to a finite 
renormalization of the T = parameters of S'meta while 
thermal fluctuations are controlled by a 'dangerous irrel- 
evant operator'. The mathematical consequence is that 
the solution of the scaling equations may be written so 
the effects of quantal fluctuations are absorbed into the 
T = parameters and only thermal effects need be ex- 
plicitly treated. We note, however, that the effects of 
quantal fluctuations are not small in general, so that a 
priori one finds the parameters (such as u) of the theory 
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are not well estimated by band theory calculations. 

The results of the calculation may be summarized as 
follows. At S = (i.e. parameters tuned so that the 
material is at the metamagnetic quantum critical point) 
as T — > the differential susceptibility dm/dh scales as 
^1/3^-2/3. ag rp _^ q t j ie S p ec jfi c ]j ea ^ coefficient 7 = C/T 

is proportional to ln/i -1 in d = 3 and to ft -1 / 3 in d = 2 
and the resistivity p(T) has the leading T dependence 
p{T) - p(T = 0) = AT 2 with A varying as ft," 1 / 3 in 
d = 3 and as hr 2 ^ in rf = 2. The crossover to the 
thermally dominated regime occurs at T ~ ft 1 / 2 (d=3) 
and T ~ ft, 2 / 3 (d=2). If 5 > then the scaling in h is cut- 
off when ft 2 / 3 ~ (5 and there is no phase transition in the 
h, t plane. If 8 < then a first order transition occurs as h 
is varied at T = 0; a line of first order transitions extends 
upwards in the h, T plane and terminates at a critical 
end-point temperature T* ~ 5 z /( d + z - 2 ) . Finally, we note 
that corrections to scaling may be numerically important, 
as will be seen from the numerical results below. 

The proceeding considerations were generic. It is 
possible to proceed further in the particular case of 
Sr3Ru207, because it seems (see below) that at ambi- 
ent field this material is very near to a weakly first or- 
der ferromagnetic-paramagnetic quantum phase transi- 
tion. The physics over a wide range of fields and temper- 
atures should therefore be describable by a generalized 
Ginzburg-Landau action for a three-component order pa- 
rameter cj> (corresponding to long wavelength fluctuations 
of the magnetization) 

f d d x f I 2 r 

So = S dyn + / —fdr < -f [V b <Pa(x, r)] + -<j> a (x, r) 

+ \ u ab(fi4>l + ^v abc 4> 2 a (j)l4> 2 c - g e gp B h ■ ${x, t) + . . . j .(5) 

Here repeated indices are summed, the ellipsis denotes 
higher order terms and the notations are as above, except 
that we have added a sixth order term and here the pa- 
rameters r, u a b, v a bc have dimension of energy. The data 
(isotropic susceptibility as T — > but some angle depen- 
dence in high-field magnetization and higher temperature 
X) require a breaking of rotational invariance in the mode 
coupling terms, but not in the quadratic one. We take <f> 
to be a dimensionless magnetization variable measured in 
units of the putative saturation magnetization 2ps /Ru 
(the important electrons are d electrons, of which there 
are four in the t2 g orbitals, leaving two holes, and the g 
factor should be close to 2). Scaling is as described previ- 
ously, except the sixth order term renormalizes the fourth 
order one and in the presence of a field the 'mass' (coef- 
ficient of the quadratic part of the fluctuations) becomes 

anisotropic, with the component corresponding to fluc- 
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tuations along the field becoming r e g = r + 3u(j) + 5v(f> . 
A Heisenberg-AF or Heisenberg-Ising crossover occurs 
when then larger of r e g or r passes through unity and 
scaling stops when the smaller of the two becomes of or- 



der unity. In the Heisenberg regime extra 'precession' 
terms in the dynamics may be important. A detailed 
analysis of the behavior of this model will be presented 
elsewhere Jf6| , extending the important work of Yamada 
and collaborators who showed that such an analysis 
was possible but did not consider the precession terms 
or anisotropic scaling and also used a simplified version 
of the 'SCR' theory instead of the renormalization group 
method. At T = 0, one may use mean-field theory as be- 
fore, provided one interprets the parameters in Eq. ^| as 
renormalized parameters. If 9u 2 /20t> > r > and u < 
(for simplicity we do not write the directional subscripts 
here) the model has a T = metamagnetic transition. 
The point r — 9u 2 /20v corresponds to the quantum crit- 
ical end-point. At the quantum-critical end-point the 
magnetization m = ((f>) and magnetic field H* are, for 
fields in the c-direction, 
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10v rn 
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Fig. I of Rcf. H shows that at low T and low ap- 
plied field the susceptibility is about 0.025^B/r imply- 
ing r w 160/^b — T w I00K. This small value implies 
a very large enhancement of the susceptibility over the 
band value, as noted previously, and implies that the 
material is near a paramagnetic-ferromagnetic transi- 
tion. For fields directed along the c axis the observed 
metamagnetic transition occurs at a magnetization of 
about 0.25 - 0.3jU B /Ru implying u cc = 3000 - 4300K 
and v ccc = 40, 000 — 80, 000K with the larger values cor- 
responding to the smaller m. The consistency of these 
estimates may be verified by substitution into Eq. ^; use 
of <7eff = 2 yields an estimate of 5 — 6T for the meta- 
magnetic field, in the range found experimentally. Ex- 
pansion of Eq. [| about the metamagnetic point yields 
Eq. with E c = -2u cc = 6000 - 8000K. The dimen- 
sionless critical field gpsH*/u ~ 0.001 so we should be 
concerned with variations which are small relative to this, 
i.e. with h ea I0~ 4 . At present rather less information 
about the spin fluctuation frequencies is available; we 
therefore normalize our results to the temperature To at 
which the differential susceptibility at the critical field is 
equal to the zero-field zero temperature susceptibility, i. e 

^(5 = 0,T=T )=x(H = 0,T = 0)=xo- 

We now present the results of a numerical solution of 
the scaling equations. Fig. || shows the h dependence of 
the differential susceptibility for several values of T, ob- 
tained in the two dimensional case using parameters rea- 
sonable for Sr3Ru207. The inset shows the temperature 
dependence of the differential susceptibility for different 
h. Note the non-monotonic temperature dependence for 
fields different from h = if the control parameter is 
tuned to criticality. Fig. || shows the specific heat coef- 
ficient 7 = C /T; in this quantity the crossover is much 
less sharp, in part because a 2d nearly critical Fermi liq- 
uid has a specific heat coefficient 7 ~ A + BT with both 
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FIG. 2: Differential susceptibility, x[j~ (c*ra/3/i), as a func- 
tion of applied field H at temperatures T/Tq = 0.05, 0.1, 0.2, 
for a two dimensional metamagnetic critical point. In- 
set: Dependence of Xo~ 1 (^ m / | 9' 1 ) on temperature T at h — 
.01, .02, .04, .08. (Normalizations discussed in text.) 
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FIG. 3: Dependence of specific heat coefficient C/T on tem- 
perature T for h = 0.01, .0.1, 0.2, 0.4 calculated for a two di- 
mensional metamagnetic critical point. Inset: Dependence of 
resistivity exponent d\n p/dlnT on T for h/H* — (lower 
curve) and 0.1 (upper curve). (To defined in text.) 

A and B divergent as the critical point is approached. 
This is an example of the corrections to scaling men- 
tioned earlier. The inset shows the resistivity exponent 
a = —d In p/d\nT plotted against temperature for h = 
and h = 0.1. The high-T resistivity exponent is not pre- 
cisely 4/3 because of the logarithmic corrections alluded 
to earlier. The crossover to the expected low-T T 2 be- 
havior is very sharp. 

In conclusion, we have presented a theory of metam- 
agnetic quantum criticality in metals, which should be 



amenable to detailed experimental tests. The universal- 
ity class was identified, a form for the order parameter 
dynamics was obtained, and detailed results were pre- 
sented for a range of physical quantities. Subsequent 
papers [jl6| will present details omitted here as well as 
results for three dimensional materials, and a comprehen- 
sive analysis of Sr 3 Ru207 and MnSi. The key assump- 
tion is that electronic degrees of freedom may be inte- 
grated out, and the fluctuating Fermi surface expressed 
in terms of an overdamped bosonic magnetization vari- 
able. Experimental tests of the theory will show whether 
this assumption, which underlies much recent work on 
quantum phase transitions in metals and indeed on non- 
Fermi-liquid physics, is justified. 
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